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ABSTRACT 



This thesis studies the treatment and evacuation of 
casualties in a Low-Intensity Conflict. The study addresses 
the effects of the evacuation policy at the medical battalion 
level, or echelon. The evacuation policy assigns a threshold 
number of days to guide surgeons in deciding which casualties 
should be treated and which should be evacuated. The main 
tradeoff in choosing this policy involves either treating and 
returning casualties to duty, or evacuating casualties to 
maintain a reserve capacity for unexpected casualty arrivals. 
The study develops simple analytical models which provide 
basic quantitative measures of the effects of the evacuation 
policy, and it establishes a framework around which more 
complex models may be built. The potential of the model is 
illustrated by examining a specific scenario. This analysis 
provides the decision maker with estimates of the number of 
available beds, the proportion of patients returned to duty, 
and the probability of being able to accommodate a mass 
arrival of casualties. 
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THESIS DISCLAIMER 

The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases of 
interest. While every effort has been made, within the time 
available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I. 



INTRODUCTION 



This thesis studies the process of casualty treatment and 
evacuation by a Naval medical battalion in a Low-Intensity 
Conflict (LIC) . The analysis focuses on the optimal policy 
for evacuating casualties. This chapter presents the 
background for understanding the military and medical 
scenarios. Following chapters fully define the problem, 
develop appropriate models, and present the results of the 
model as applied to a specific decision-making scenario. 

A. PROBLEM SETTING 

In the late 1980's and early 1990's the Marine 
Expeditionary Force (MEF) has been deployed as the initial 
force in United Nation (UN) peace-keeping and relief efforts. 
The MEF, consisting of a Marine division and its support 
services, may deploy to foreign and hostile environments for 
possibly two to three months. Support services for the MEF 
include a Naval medical battalion to provide medical support 
for approximately 15,000 marines. 

1. Military Scenario 

The type of combat environment expected for this type 
of operation can be termed a Low-Intensity Conflict (LIC) . In 
the UN peacekeeping role, the MEF may deploy to an area where 
the greatest medical threat is that of local disease. 
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Marines might also expect a low rate of combat injuries. 
Beyond these expected disease and battle casualties, there is 
typically the potential for an unexpected episode resulting in 
many casualties. 

2. Medical Conditions , 

All casualties can be placed into two distinct 
classifications: Wounded in Action (WIA) , and Disease and 
NonBattle Injuries (DNBI ) . For medical purposes, the arrival 
and treatment rates for these casualty classes are considered 
separately. In the LIC environment, arrival rates of DNBI and 
WIA casualties should normally be low. However, a large 
casualty occurrence could result in sudden, but short-lived, 
increase in the WIA rate. 

In this scenario, the medical planners' decisions 
include the determination of an evacuation policy. The 
evacuation policy stipulates the maximum amount of time that 
a patient may remain at a treatment facility before he should 
be evacuated. This policy implicitly assumes that if patients 
are allowed to remain in beds for a longer period of time, 
then more of these patients will be returned to duty. 
However, the resulting larger population of recovering 
patients reduces the facility's ability to accommodate a 
sudden influx of wounded casualties. 

Under the conditions of this scenario, medical 
planners should set an evacuation policy that allows for the 
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accommodation of an unexpected and large casualty "spike." 
Simultaneously, planners should implement a policy that 
results in the return to duty of most patients during day-to- 
day operations. This fundamental tradeoff forms the basis of 
the evacuation policy problem. 

B. THE EVACUATION POLICY PROBLEM 

The goal of this thesis is to model the casualty treatment 
and evacuation process, and to study quantitative methods to 
aid decision makers in setting the evacuation policy. 

The evacuation policy sets a threshold number of days for 
patient treatment . If the evaluation of an incoming 
casualty indicates a treatment time less than this 
threshold, then the casualty is treated in anticipation of 
return to duty. Otherwise, the casualty is stabilized and 
evacuated. Henceforth, this treatment threshold is 
synonymous with the evacuation policy. 

Analysis is limited to a LIC scenario with normally low 
casualty rates that are within the capacity of the medical 
battalion's resources. Based on this constraint, we can 
analyze the impact of the evacuation policy on the medical 
battalion's effectiveness. 

This thesis begins by discussing the problem's parameters 
and tradeoffs, and identifying appropriate assumptions and 
measures of effectiveness. The study continues by developing 
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the framework for a real time decision aid based on analytical 
techniques. The resulting model is used to analyze a decision 
problem supplied by the Naval Health Research Center (NHRC) to 
illustrate the potential usefulness of these analytical 
techniques. Finally, the study addresses some modeling 
shortcomings and possible alternate models. 
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II. PROBLEM DISCUSSION 



A. ASSUMPTIONS AND PARAMETERS 

In approaching this problem, the first step is to define 
the system of interest, namely the medical battalion. 
Figure 1 is a flow diagram of the general casualty treatment 
and evacuation process. The medical battalion's position in 
the medical evacuation chain is evident in this figure. 




Figure 1 General Casualty Treatment Flow Diagram. 



To concentrate analysis on the medical battalion, several 
assumptions are made to reduce Figure 1 to an uncomplicated 
representation of the medical battalion. The resulting 
structure consists of the medical battalion, its input of 
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casualties, and its output of patient evacuations and returns 
to duty. 

From this point on, the term casualties is used to 
describe those wounded or injured personnel that have not yet 
received treatment by the medical battalion. Patients , on the 
other hand, are those personnel that have been admitted to the 
medical battalion for treatment. 

1. Casualty Arrivals 

As indicated in Figure 1, the medical battalion 
receives casualties direatly from the' field or by way of 
battalion aid stations. Experience has shown that emergency 
treatment by company medics or battalion aid stations greatly 
increases a casualty's chance of surviving and returning to 
duty [Ref. 1]. This study assumes that this preliminary 
treatment is reflected in the model for the distribution of 
casualty treatment times. This assumption allows us to 
disregard the point of origin of arriving casualties. 

Besides the point of origin, consideration should be 
given to the types of arriving casualties. Ideally, 
casualties can be categorized in terms of specific types of 
wounds, diseases and injuries. For some specific scenarios, 
in which a limited number of ailments are expected, this 
approach might be feasible. However, this study considers 
only the military's two general casualty types: Disease and 
Nonbattle Injury { DNBI ) , and Wounded in Action (WIA) . These 
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categories are well understood by military and medical 
planners, and current procedures are in place to estimate 
casualty rates for each category. 

2 . Patient Departures 

This analysis assumes that patients leave the medical 
battalion for one of two reasons: return to duty, or 
evacuation. The focus is maintained on the tradeoff between 
casualty evacuation and return to duty by not analyzing for 
mortality effects. If desired, the small number of patient 
deaths could be addressed, possibly by including these figures 
with patient evacuations. However, mortality rates in the 
medical battalion should be low enough in this type of 
conflict so that the results are not affected. 

3. Medical Unit Operations 

A significant assumption for this study pertains to 
the clinical aspects of casualty treatment. At many points in 
the treatment of casualties, doctors make clinical 
evaluations, or forecasts, that affect the treatment process. 
Examples of this include the triage of incoming casualties and 
the day-to-day evaluation of patient conditions. 

Triage is "... the sorting and assignment of treatment 
priorities to various categories of wounded." [Ref. 2] The 
triage decisions determine whether or not a casualty is 
admitted for treatment, stabilized for evacuation, or 
immediately returned to duty. Thus, the triage officer is 
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forecasting the required treatment time and effort required 
for each casualty. For this study, triage and all other 
clinical forecasts are assumed to be made perfectly. 

From a logistics standpoint, this analysis does not 
account for limited transportation or supply resources for the 
medical battalion. The model assumes that assets are always 
available to transport casualties, and supply support is 
always sufficient to support full bed utilization. The model 
leaves it to the medical commander to decide the required 
amount of transportation and supply support based on the 
calculated measures of effectiveness. 

B. POLICY TRADEOFFS 

This study assumes that the user of the model is an expert 
decision maker. In other words, the decision maker 
understands the tradeoffs that must be made in providing the 
best medical service for a Marine division. A discussion of 
these tradeoffs underscores the importance of the evacuation 
policy. 

1. Mission 

An understanding of the medical battalion's mission 
is required to grasp the significance of these tradeoffs. 
This mission, according to Emergency War Surgery, follows: 

The ultimate goal of combat medicine is the return of the 
greatest possible number of soldiers to combat and the 
preservation of life and limb in those who cannot be 
returned. [Ref 2] 
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Clearly the evacuation policy greatly affects this 
mission. This policy directly influences the medical 
battalion's ability to fulfill both aspects of its mission. 

2. Primary Effects 

By setting the evacuation policy at a large value, 
many patients can eventually be returned to duty. This policy 
causes more of the medical battalion's beds to be occupied. 
However, this increased bed utilization diminishes the ability 
of the medical battalion to respond to a large surge in 
casualty arrivals. 

A large group of casualties arriving to find a full 
medical facility results in casualties either waiting for 
treatment, or being evacuated unnecessarily. This situation 
is undesirable for the following reasons: 



. . . those casualties most likely to return to duty will be 
sent farthest away, since they are the most transportable 
and in the best condition to tolerate evacuation. Those 
casualties least likely to return to duty will be kept in 
forward areas, due to their serious and unstable 
condition. These may then pose a much greater logistical 
demand upon the system due to the serious nature of their 
wounds requiring intensive care. [Ref. 3] 



If planners choose a shorter evacuation policy, the 
resulting bed utilization is much lower. This has the 
advantage of allowing the accommodation of a large group of 
casualties. However, this policy results in fewer patients 
returning to duty. 
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Thus, the medical battalion's patient returns to duty, 
and its ability to accommodate large casualty arrivals, are 
directly affected by the choice of evacuation policy. 

3. Secondary Effects 

Besides these considerations, other tradeoffs 
accompanying the evacuation policy lead to secondary effects. 
These effects are termed secondary because of their dependence 
on the previously discussed factors. They are not necessarily 
secondary in importance to the military commander. Relevant 
factors are typically hard to measure, thus making 

quantitative analysis of their effects difficult. Some of 
these factors include: 

• Higher bed utilization requires increased medical supply 
support . 

• Increased evacuation rates result in greater 

transportation requirements. 

• Increased evacuation rates also generate an inflow of 
inexperienced replacements to the combat units. 

• Any inability of the medical battalion to handle mass 
casualties results in unnecessary casualty evacuations and 
treatment delays . 

• This inability also reduces the confidence and morale of 
the combat units. [Ref. 1] 

These are only a few of the many effects that a 
medical planner may consider in choosing an evacuation policy. 
The importance of each of these factors varies depending on 
the actual military situation. Consequently, we define a few 



10 



important measures that may help the planner estimate the 
actual effects of an evacuation policy. 

C. MEASURES OF EFFECTIVENESS 

Many concepts used in queueing theory are applicable to 
the evacuation policy. This discussion presents some of these 
concepts as they pertain to the evacuation policy problem. 

1 . Primary Measures 

The medical battalion's mission effectiveness is 
primarily measured by the flow of patients returning to duty 
and the ability to treat the majority of arriving casualties. 
In evaluating the medical battalion's ability to treat a large 
influx of casualties, it is assumed that these mass casualties 
arrive as a group 1 . Thus, the maintenance of a reserve bed 
capacity becomes important to the planners, 
a. Patient Returns to Duty 

The first measure of patient returns to duty is the 
average, or mean, rate of return to duty, X RTD . This measure 
provides a steady-state approximation to the average number of 
patients returned to duty each day. This measure can assist 
the medical planner in estimating requirements such as 
transportation assets. 



This assumption is based on a large casualty influx that 
lasts for a short time (i.e., < 24 Hours) as compared to the 

campaign time (i.e., « 60 to 90 days). 
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Alternatively, the decision maker may be more 
concerned with the expected proportion of patients returning 
to duty, p ?Tr . This measure provides a more direct indication 
of the effectiveness of the medical battalion in returning 

b. Resen r e Capacity 

The most straightforward estimate of the medical 
battalion's reserve capacity is the expected, or average, 
number of available beds. This provides the medical planner 
with a measure of long-term effects such as required supply 

Another useful measure is the probability of having 
at least a specified number of beds available. This provides 
the planner with an actual probability of accommodating a 
large casualty surge. 

2, Secondary Measures 

Beyond these primary measures, some secondary measures 
cf effectiveness further illustrate the effects of the 
evacuation policy. Some of these measures include: 

• The expected rate and number of patient evacuations, 

• The expected rate and number of replacements required, and 

• The expected number of beds occupied. 

By using analytical models, many of these primary and 
secondary measures of effectiveness can be calculated using a 
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hand-held programmable calculator. Alternatively, the results 
can be plotted, for a range of evacuation policies, to provide 
a graphical representation of the measures of effectiveness. 



In either case, this analysis results 
quantitative measures to assist the medical 
with the assignment of an evacuation policy. 



in accessible 
decision maker 
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III. MODEL STRUCTURE 



We now develop an analytical model for the treatment and 
evacuation process. An analytical model, which captures the 
essential components of the real problem, offers the advantage 
of precise, easily calculated measures given some reasonable 
modeling assumptions. 

Figure 2 is a block diagram of the treatment and 
evacuation model. Development begins with a discussion of the 
flow parameters for the model, and then addresses the patient 
treatment time considerations. 



Patient r — N 

Characteristics 1 — 1/ 



Casualty i — N 
Rates — v 



Patient 
Row Paths 






Casualty 

Flow 

Block 



Treatment 

Parameters 



Admission 

Parameters 



Treatment 

Time 

Block 



[Evacuation 
and 
Return to 
Duty 
Measures 



Bed 
Utilization 
and 
Reserve 
Capacity 



Figure 2 Model Block Diagram. 
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It is useful to define several terms and indices that are 
used consistently in the following model development. Some of 
these include: 

Casualty Types (Referenced with the subscript i) : 

• W = WIA = Wounded in action, and 

• D = DNBI = Disease and nonbattle injuries. 

Treatment Types (Referenced with the superscript j): 

• A = Admissions treated for return to duty, and 

• E = Evacuees stabilized before evacuation. 

Other (Superscript) : 

• R = Casualties immediately returned to duty. 



A . CASUALTY FLOWS 

The flow parameters for the casualty treatment model 
include patient characteristics, casualty rates, and patient 
flow paths. These factors affect the medical treatment system 
and consequently influences our measures of effectiveness. 

1. Patient Characteristics 

The most important patient characteristics concern the 
type and severity of arriving casualties' wounds or injuries. 
Accurate identification and triage of the various casualty 
severities are crucial to medical operations. 

During the triage of incoming casualties, medical 
commanders can readily figure out the proportion of casualties 
who are either immediately scheduled for evacuation. 
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immediately returned to duty, or admitted to the hospital for 
treatment in anticipation of return to duty. As mentioned 
earlier, these forecasts by the triage officer are assumed to 
be made perfectly. 

By placing arriving^ casualties in one of these three 
groups, we can define the following proportions for i = w and 
D (DNBI and WIA casualty types) 

pf = proportion of arriving type i casualties scheduled 
for evacuation, 

Pi = proportion of arriving type i casualties immediately 
returned to duty, and 

Pi = proportion of arriving type i casualties admitted to 
the medical battalion for treatment; 

where pf + pf + pf = 1 . 

2 . Casualty Rates 

Arrival rates for DNBI and WIA casualties are also 
essential inputs to the model. Casualty rate estimates are 
normally calculated during the planning process of a military 
operation. These casualty estimates should be used by 
commanders to make initial evacuation decisions. However, 
once the operation begins, actual casualty rate data should be 
used as a basis for updating the evacuation policy. 

Casualty rates are normally expressed as the number of 
casualties arriving per one thousand combat troops per day. 
For this analysis, we convert these casualty rates into the 
following daily rates: 
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X D = the average number of DNBI casualties arriving per 
day , and 

X w = the average number of WIA casualties arriving per 
day . 

Note that these are estimated rates of casualties 
arriving at the medical battalion for triage. These are not 
necessarily the rates of patient admission. Patients admitted 
to the medical battalion include those being treated for 
return to duty, and those being treated before evacuation. 

3. Patient Flow Paths 

At this point, it is worthwhile to identify further 
the different flow paths in the patient treatment model. 
Figure 3 illustrates these paths as they are determined during 
the triage process. Notice that some casualties are returned 
to duty immediately during this screening, and the remaining 
casualties are admitted to the hospital. 

As Figure 3 shows, there are four patient streams 
being admitted for hospitalization. The DNBI and WIA 
casualties are shown as distinct flows paths in the triage 
process. Each of these paths is further split into those 
patients being treated for return to duty and those being 
treated before evacuation. 
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Casualties are assumed to arrive according to a 
Poisson process. Thus, by decomposition, the casualty 
admission rates for each of these paths can be found by 
multiplying the casualty arrival rates by the patient 
proportions : 

=p/*i, V(i,j), (l) 

where is the admission rate of casualty type i for 

treatment type j . 
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B. PATIENT TREATMENT TIMES 



We now consider the patient treatment times for each of 
the patient flow paths. This discussion addresses the 
decisions affecting the evacuation of arriving casualties, 
then proposes a model for the patient treatment times. 

1. Evacuation Decisions 

There are two decisions affecting the evacuation of 
casualties. These are the evacuation policy and the 
evacuation schedule. The evacuation policy, e p/ was 
previously defined as the maximum amount of time that a 
casualty may be treated before he should be evacuated. For 
those patients identified as evacuees, the evacuation 
schedule, e s , is an estimate of the amount of hospitalization 
time needed before their evacuation. 

The evacuation schedule accounts for the time needed 
for some patients to stabilize following emergency surgery. 
The actual treatment time, for these evacuees, is a random 
variable. However, previous analysis has shown that the 
evacuation schedule can be approximated by a fixed value 
[Ref. 4]. Separate evacuation schedules may be set for DNBI 
and WIA patients; however, this analysis assumes a single, 
fixed evacuation schedule for all evacuating casualties. This 
value represents expected delay time before evacuation for 
these casualties. Historical estimates provide the evacuation 
schedule for this analysis. 
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The evacuation policy is the focus of this study. It 
is assumed to apply equally for DNBI and WIA casualties. The 
qualitative effects of the evacuation policy have been briefly 
discussed. Succeeding analysis concentrates on the 
quantitative effects of this policy on the medical battalion's 
effectiveness . 

2 . Treatment Time Model 

For the purposes of this study, treatment time is the 
number of days of treatment needed before allowing a patient 
to return to duty. For the treatment time model, we consider 
only those casualties admitted to the medical battalion in 
anticipation of return to duty. This allows us to more 
accurately depict the bed utilization of the medical 
battalion. A representative data set, composed of treatment 
days required to return patients to duty, forms the basis for 
the distribution model. These data, provided by the Naval 
Health Research Center, represent typical casualty treatment 
times for a LIC scenario. Table I presents these data for 
both DNBI and WIA casualties. 

Table I data indicates the number of instances of 
treatment days required for return to duty of casualties. 
This frequency distribution can be modeled with a parametric 
distribution to provide general analytical results for this 
study. One must exercise caution when using current data, as 
it was most likely collected when an evacuation policy was 
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TABLE I RAW CASUALTY TREATMENT TIME DATA [Ref. 4] 



Days 


Number 
of DNBI 
Cas . 


Number 
of WIA 
Cas . 


Days 


Number 
of DNBI 
Cas . 


Number 
of WIA 
Cas . 


: d 
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40 
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3 
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10 
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40 
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11 
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10 
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11 


42 
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13 


53 


9 


43 


1 


0 


14 


57 


11 


44 
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0 


15 


43 


8 


45 
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0 


16 


37 
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40 
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17 


22 


10 
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14 


27 
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46 
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19 
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22 


11 


50 


1 


1 


21 


16 


7 


51 


0 


0 


22 


20 


t 


52 


1 


0 


23 


20 


8 


53 


c 


1 


24 


20 


2 


54 
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1 


25 


11 


3 


55 


0 


0 


26 


11 


* 


56 


o 

u 


0 


27 


7 


* 


57 


1 


0 


28 


8 


2 


58 


0 


0 


29 


4 


2 


| > 59 


6 


608 




| Totals 


1930 


1616 
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being used. Thus, the upper tail of the distribution tends to 
be censored. Despite this property, the data provide a 
starting point for comprehending the form of the distribution 
of casualty treatment times. 

a. Data Discussion , 

The treatment data set not only provides an 
indication of the distribution of treatment times, but it also 
provides an estimate of the patient characteristics. The data 
includes instances of patients who, requiring zero days of 
treatment, may be immediately returned to duty. In addition, 
it includes patients requiring more days of treatment than the 
maximum permitted evacuation policy of 60 days. Analysis of 
these outliers allows us to estimate the proportion of 
patients immediately returned to duty. 

Those patients requiring zero days of treatment 
typically suffer from minor injuries or battle fatigue and do 
not require hospitalization. Although these arrivals cause 
short-term stress on the system, they do not affect bed 
utilization. These casualties may be treated, rest in holding 
wards for a few hours, and then are returned to duty. 

The data in Table I suggests that the following 
proportions of arriving casualties not requiring 
hospitalization before return to duty: 
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* _ 97 

" TciT 



1930 



= 0 . 050 , 



( 2 ) 




= 0.402 . 



These figures may be compared to those recommended by Marine 
Corps planners of 20% for DNBI and 50% for WIA casualties 
[Ref. 6] . In reality they are heavily dependent on the tempo 
and nature of operations, and may vary over a wide range. 
Fortunately, medical battalion personnel can easily monitor 
these figures and apply the measured values to actual 
situations. These data elements, representing minimally- 
injured casualties not admitted to the hospital, provide an 
estimate of the proportion of patients that will be returned 
to duty immediately. 



represents casualties with severe injuries and with little 
hope of return to duty. This group is likely to be a mixture 
of casualties with stable, but debilitating, injuries and 
severe life threatening injuries. In either case, the medical 
battalion should not treat these casualties intent upon their 
return to duty. Since the maximum allowed evacuation policy 
is 60 days, only those casualties with treatment times less 
than 60 days are considered to have any potential for return 
to duty. These casualties are subject to the evacuation 
schedule; and, will burden the system until their condition 
permits evacuation. 



The opposite end of our treatment data spectrum 
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These parameters can be expressed in terms of the 
previously defined casualty proportions, or equivalently in 
terms of the cumulative distribution function for treatment 
times. Figure 4 illustrates this point with a plot of the 
cumulative treatment times -for the DNBI data of Table I using 
a typical evacuation policy of 20 days as an example. 




Figure 4 DNBI Cumulative Treatment Time Data. 



In this figure, F D (t) is the cumulative distribution function 
of treatment time for DNBI casualties. The plot clearly shows 
the following relationship between the triage proportions and 
F D ( t ) given an evacuation policy of 2 0 days: 
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Pd = 1 " F D (20) , 

Pd = F d { 20) - F d ( 0) , and 
Pd = ^(0) . 



( 3 ) 



The relationships illustrated in Figure 4 and tabulated in 
Equation (3) also hold for the WIA casualty type. 

b. Distribution Fitting 

The remaining data in Table I represents the 
treatment times of casualties that may be admitted to the 
medical battalion with reasonable expectation of their return 
to duty. These data are used to estimate an approximate 
distribution of treatment times for all return-to-duty 
patients. WIA and DNBI data from Table I were analyzed 
separately using the distribution fitting tools in AGSS 2 . In 
both cases, a lognormal distribution provides a better fit 
than alternatives such as the exponential and Weibull. 
Figure 5 shows the DNBI data from Table I for treatment days 
ranging from 1 to 59 . This figure illustrates the lognormal's 
more accurate approximation of the distribution of treatment 
time data as compared to the other two distributions. 

For the DNBI data, the resulting lognormal 
distribution had a chi-square goodness of fit significance 



A Graphical Statistical System (AGSS) is IBM's commercial 
equivalent to the beta test package, GRAFSTAT, used at the Naval 
Postgraduate School . 
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Probability Distribution Fits 
DNBI Treatment Times 




Figure 5 DNBI Data Histogram and Distribution Fits. 

level of 0.042. The fit of the WIA data resulted in a 
significance level of 0.0054. Appendix A contains the details 
of this data analysis. The resulting lognormal distributions 
are accepted as reasonable approximations of the WIA and DNBI 
data in the range from 1 to 59 treatment days. 

C . MODEL SUMMARY 

At this point, the treatment and evacuation process has 
been dissected into its basic patient flow and treatment 
components. Figure 3 provides a good summary of these 
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components. The figure clearly shows the four patient flow 
paths as they are admitted to the medical battalion for 
treatment. The proportion of patients along each path is 
determined by the distribution of the casualty severities and 
the evacuation policy. 

The distribution of treatment times, resulting from the 
analysis of Table I data, is the mixed distribution, defined 
by: 

e 

F 1 (t) = F 1 (0) +F r 1 (0)Jg 1 (c)dt,^i l ( 4 ) 

0 

where g^t) is the lognormal probability distribution fitted 
in the previous section. 

These modeling elements may now be blended to determine 
several measures of effectiveness for the medical battalion. 
The resulting measures aid the decision maker in selecting the 
appropriate evacuation policy for the medical battalion. 
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IV. MODEL OUTPUTS 



A. BED UTILIZATION 

Several measures of bed utilization are developed in this 
section by combining the treatment time model with the input 
flow parameters. The average, or mean, number of available 
beds is the initial measure of bed availability. Often, an 
estimated probability of accommodating a group of casualty 
arrivals may be more useful than a simple expected-value 
estimate for available beds. For this measure, we must 

approximate the probability distribution of bed occupancy. An 
Erlang-B loss model is developed for this purpose. Finally, 
the characteristics of the Low-Intensity Conflict allow the 
use of a simpler Gaussian model to approximate the 
distribution of bed occupancy. 

1. Expected Value Approach 

Assuming a maximum bed capacity, k, for the medical 
battalion, the expected number of available, or empty, beds is 
the maximum bed capacity minus the expected number of beds 
occupied, or 

E[N av ] = k - e[iv oc ] (5) 
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where E[N oc ] is the expected number of occupied beds, and 
E [N av ] is the expected number of available beds. Thus, we need 
to estimate the number of occupied beds, E[N oc ] . 

As the flows paths in Figure 3 show, bed occupancy is 
made up of patients from four separate flow paths. 
Considering the contribution to bed utilization from each of 
the four paths, E[N oc ] can be found by 

E[N oc ] = £ E[N±] = E[Nd ]+E[N$] +E[n£] +E[N%] , ( 6 ) 
VU'j) 

where E [Nj] is the expected number of beds occupied by type i 
casualties undergoing type j treatment. 

Each term in Equation (6) can be found by applying 
Little's result. This result 

. . . states that the average number of customers in a 
queueing system is equal to the average arrival rate of 
customers to that system, times the average time spent in 
that system. [Ref. 7] 

For the medical battalion, this implies that the average 
number of occupied beds is equal to the average admission rate 
of patients, multiplied by the average treatment time for 
those patients. 

In applying Little's result, this analysis assumes 
that the medical battalion is never full and casualties do not 
queue. Thus, average time in the system is the same as the 
average treatment time. With this assumption, applying 
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Little's result to each term in Equation (6) yields the 
following equation: 

= £ (A^Bls/] ) , (7) 

vu.j) 

where X{ is the average arrival rate, and E[Sj] is the average 
treatment time, for each casualty type i and treatment type j . 

The terms in Equation (7) can be found based on the 
parameters already introduced. Given the casualty 

proportions, estimates of the arrival rates may be calculated 
from the casualty arrival rates as suggested by Equation (1) . 
Furthermore, the expected treatment time for evacuating 
patients depends only on the evacuation schedule, which is 
assumed constant for this analysis: 

S[S C B ] = Elsf] = e s . (8) 

This leaves E[Sj] and E[Sj] as the only still unknown terms in 
Equation (7) . 

Recalling the casualty treatment time distributions, 
from Equation (4), defined in terms of g^t), we can find the 
expected treatment time for patients being treated in 
anticipation of returned to duty by: 

BIS/] =Bj( 0) f f *\ lt) dt, Vi. (9) 

G i (e V 
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Considering this relationship, the terms in 
Equation (7), for return-to-duty patients, can be written as: 



A *-E[S?) = \ 1 -Pi-F i (0) ft 
0 * 



sr.it:) 

G ii e p ) 



dt, Vi. 



( 10 ) 



But, by recognizing that 



G 1 (e p ) = p? = F d (e p ) - F 1 (0) , Vi, 
Equation (10) can be reduced to: 

k.-E[S?] = X.-FliO) I fg.it) dt, Vi. 



( 11 ) 



( 12 ) 



We now have all the terms needed to calculate the 
expected number of available beds in the medical battalion. 
The expected number of occupied beds becomes 3 : 



£[««;] = W°> / t-g D (t)dt 
0 * 

e P 

+ *.„-F'„( 0 ) f t-g w (t)dt 
0 * 

+ U| + *£> 



( 13 ) 



where the flow rates X D and are given, and the evacuation 
decisions e p and e s are constants. The expected number of 



This result assumes that DNBI and WIA casualties are 
sorted based on the same evacuation policy. If desired, separate 
evacuation policies could certainly be used. 
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available beds is found by subtracting the result from 
Equation (13) from the total bed capacity of the medical 
battalion (See Equation (5) ) . 

This measure of the expected number of available beds 
may not be sufficient for the medical decision maker. Often, 
he may wish to know the probability of at least a specific 
number of beds being available. This type of measure requires 
a more detailed treatment model. 

2 . Erlang-B Loss Model 

The characteristics of the medical facility lead to 
the Erlang-B loss model. This model, originally for telephone 
switching systems, provides an initial off-the-shelf but 
reasonably appropriate analytical model for the present 
situation . 

The Erlang-B model (Also termed the Erlang loss 
system) is equivalent to the M/M/k/k service system, in which 
the M/M/k/k notation describes a service facility receiving 
Poisson arrivals, with exponential service times, and k 
servers with no holding or queueing capacity. 

In this system, customers (i.e., casualties) arrive 
according to a time -homogeneous Poisson process. If at least 
one of the k servers (i.e., hospital beds) is available then 
the customer enters the system and spends an exponentially 
distributed time in service. If all of the k servers are busy 
(i.e., all beds are full), the system is at capacity. A 



32 



customer who arrives when all of the servers (beds) are 
occupied is lost from the system. 

For this model, the steady-state probability of being 
in any state n, namely P n , can be found by evaluating the 
formula : 



Pn 



{X/\j) n /n\ 

E U/|i)Vii 



0 , 1 , 2 , 



• ,k , 



( 14 ) 



where 

X is the rate of the Poisson arrival process; 

|i is the rate of customer (casualty) service; 
n is state of the system; and 

k is the number of servers (beds) in the system. 



If we define S as the service time, then the expected 
or average service time, E[S], is 1/|X. By substituting for [I 
in Equation (14) it follows that 



Pn = 



(Vg[S] ) n /n! 

£ U\E[S] ) i /i\ 

i*0 



rx - 0,1,2, 



. ,k. 



( 15 ) 



The following is an important fact: if we now consider 
the system M/G/k/k, where the service times have any general 
distribution (denoted by G) , it can be shown that 
Equation (15) continues to hold for this more general system. 
It is said that the Erlang-B formula is insensitive to the 
form of the service time distribution, being dependent only on 
the mean service time. The resulting model is the Erlang loss 
system. [Ref. 8] 
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To apply these results to the medical battalion, we 



must first define the combined arrival rate, X, T : 

k T ~ ^ D + 



and the combined mean service time, E[S T ]: 



E[S t ] 



ElN t pg] 



( 16 ) 



( 17 ) 



where E[N oc J is found by Equation (13). Substituting these 
results into Equation (15), we find that: 



PrT P^ 0 c=n) 



(X r -g[S r ])Vn! 

£ (A. T -£tS T ] ) i /i ! 
i=o 



0 , 1 , 2 , 



• ,k. ( 18 ) 



Since N oc = k - N av , the probability of at least n beds being 
empty, P(N av > n) , is derived by: 

P{N av zn) = P(k-N oc zn) 

k-n-l ( 19 ) 

= P{N oc <k-n) = £ P d . 

1=0 

The insensitivity of this result to the treatment time 
distribution makes the Erlang-B model a robust model for 
analyzing the bed occupancy level of the Medical Battalion. 

3 . Normal Approximation 

The nature of this problem allows us to apply a normal 
approximation to the Erlang-B model results. In the LIC 



34 



scenario, we are assuming that the steady-state arrival rates 
are within the capacity of the medical battalion. This 
assumption allows the medical battalion bed system to be 
modeled as an infinite server queue. 

This assumption implies that the distribution of p n in 
Equation (18) is approximately Poisson with parameter X T *E[S T ] , 
or E [ N oc ] . Further, since the expected number of beds occupied 
is greater than 20 for the scenario of interest, the 
distribution for p n can be approximated as normal. 

Thus, given the expected value of the number of 
occupied beds, E[N oc ], the distribution for N oc is 
approximately normal with mean and variance equal to E[N oc ]. 
Using this fact, we can find the P(N av > n) in terms of the 
standard normal random variable z as follows 



P(N av zn) = P ( k-N oc z n) 

= P(N oc z k-n) 



= P 



Kc-E[N oc ] ^ k-n-E[N 0l 



k-n-E[N ] 



L ) 



P\z <. 



f £ 



-n-E[N 0 J 






( 20 ) 



The accuracy of this normal approximation depends on 
the patient admission rates and their mean service times. 
Analysis found that the normal approximation is accurate up to 
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a bed utilization of approximately 0.90, where bed 
utilization , p, is 

p = = EtN oc ] (21) 

For the LIC, the bed utilization is typically large enough to 
allow the use of the normal approximation. At the same time, 
the bed utilization small enough to justify the use of the 
Poisson approximation to the Erlang model. 

The normal approximation provides a second and, when 
applicable, convenient method to evaluate the distribution of 
bed occupancy. Figure 6 shows a comparison between the 
Erlang-B results and those found with the normal 
approximation 4 . It is clear from this figure that the normal 
approximation is almost indistinguishable from the Erlang 
model using reasonable parameters for casualty rates, 
treatment times and evacuation policies. 

B. RETURN TO DUTY 

The other primary measure of effectiveness concerns 
patient returns to duty. Assuming perfect forecasting during 
the triage process allows a direct calculation for the patient 
returns to duty. 



The Erlang-B distribution was solved analytically using 
the algorithm found in Appendix B. 
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Number of Beds Occupied 



Figure 6 Normal Approximation to the Erlang Model. 



For a given evacuation policy, e p/ the proportion of 
patients returning to duty following treatment is the 
cumulative distribution of treatment times, for DNBI and WIA 
casualties, evaluated at e p . Assuming perfect forecasting, 
all treated patients eventually return to duty. Thus, the 
total proportion of arriving casualties returning to duty, 
PrTD' i- S : 



Prw = Pi + Pi 

= 1 - pf, Vi . 



( 22 ) 



In a more complex treatment process model, the uncertainty of 
the triage process leads to some admitted patients not 
returning to duty. 
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C. SUMMARY 

The complete casualty treatment model provides the 
decision maker with two general measures of effectiveness: 
patient return-to-duty measures and bed utilization measures. 
Each of these general results may be evaluated by medical 
planners, within the current mission and situation, to assist 
in their decision making. 

An advantage of this model is its design allowing for a 
variety of specific measures. For this study, however, the 
model results are demonstrated based on the previously defined 
measures of effectiveness. 
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V. 



MODEL RESULTS 



The models developed in the preceding chapters may now be 
applied to a specific problem supplied by the Naval Health 
Research Center. Numerical and graphical results are 
presented to illustrate the potential for these types of 
analytical models. 

All calculations and plotting for this analysis were 
performed in Microsoft Excel 4.0. A sampling of the Excel 
worksheets is contained in Appendix C. 

A. NUMERIC INPUTS 

The proposed problem considers a medical battalion with a 
540-bed capacity. This facility is providing medical support 
for a combat force of 15,000 troops. The estimated casualty 
rates for this operation are: 

DNBI rate =3.5 per 1000 strength per day, and 
WIA rate = 1.5 per 1000 strength per day. 

Beyond these steady-state rates, the enemy is believed to have 
the capability to inflict a one-time only WIA rate of 15 per 
1000 strength for one day. 

These casualty rates are converted to this study's 
notation by: 
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X D = 3.5 x 15 = 52.5 Casualties/Day, 

= 1.5 x 15 = 22.5 Casualties/Day, 

X Max = 15 x 15 = 225 Casualties/Day. 

Further, the maximum one-day increase in the number of 
patient admissions is calculated as 

Casualty Spike = - X w ) x (1 - p£) 

= {225 - 22.5) x 0.598 
= 121 patients. 

Additional numeric inputs include: 

k = 540 beds (From the problem statement) , and 
e s = 4 days (Historical estimates [Ref. 4]). 

Parameters for the lognormal models are (See Appendix A) : 
j 1 D = 2.1 and CJ D = 0.737, 

|l w = 2.25 and U w = 0.75. 

These values comprise the numeric inputs to the models 
developed in the previous chapters. Using these figures, 
baseline results are calculated for the measures of 
effectiveness. In addition, some parameters are varied to 
provide some sensitivity analysis for the results. 

B. BED UTILIZATION RESULTS 

We begin by estimating the expected number of available 
beds, E[N av ], for a range of evacuation policies. Numerical 
evaluation of Equations (5) and (13) allows us to estimate 
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these expected values. Figure 7 plots the expected number of 
beds available versus the evacuation policy. For this 
scenario, this figure suggests an evacuation policy of 14 days 
to accommodate the maximum one-day casualty arrival event, or 
spike . 




Figure 7 Results for Expected Number of Beds. 



We may now estimate the medical battalions' ability to 
accommodate the large casualty arrivals in terms of a 
confidence level. The normal approximation model was used to 
find a 95% confidence level of available beds. These results 
are shown in Figure 8. Notice in this figure that the 
indicated evacuation policy has been reduced to between 12 and 
13 days. Thus, the recommended evacuation policy would be 12 
days . 



41 




Figure 8 Results for N av with 95% Confidence Level. 



Considering these results. Figure 9 was created to examine 
the range of feasible evacuation policies for our estimated 
casualty spike. For a casualty spike of 125, which 
approximates the parameter for our problem, the evacuation 
policy should be less than 17 days. Between 17 and 10 days, 
the choice depends on the confidence level desired by the 
decision maker. For evacuation policies less than 10 days, 
accommodation of a casualty arrival of size of 125 is assured. 

These results provide recommended evacuation policies, for 
this specific scenario, based on bed utilization and a desired 
confidence level. The decision maker should also consider the 
effect that the choice of evacuation policy will have on 
casualty returns to duty. 
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Evacuation Policy (Days) 



Figure 9 Sensitivity Analysis of Confidence Level. 

C. RETURN-TO-DUTY RESULTS 

We first consider the proportion of potential patients 
that are returned to duty. This excludes those casualties 
that do not require hospitalization before return to duty. 
This exclusion gives a direct indication of the association 
between returns to duty and bed usage. 

Figure 10 illustrates the results for the proportion of 
patients returned to duty. Recall that this measure is based 
on the cumulative treatment time distribution, hence its 
monotonically increasing shape. Evacuation policies, based on 
the bed utilization results, ranged between 10 and 17 days. 
These evacuation policies lead to return-to-duty rates of 0.58 
and 0.82 respectively. The 12 -day evacuation policy, 
recommended for a 95% confidence of accommodating the casualty 
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spike, results in a return-to-duty proportion of 0.68, as seen 
in Figure 10. 




Evacuation Policy (Days) 

Figure 10 Proportion of Patients Returned to Duty. 

D . UTILITY OF RESULTS 

For this specific problem, the recommended evacuation 
policies range from 10 to 17 days. The medical planner must 
consider the military situation in selecting the threshold 
number of days for the evacuation policy. 

1. Maximizing Return to Duty 

To maximize returns to duty, the planner would like to 
set the evacuation policy at 17 days. This would allow the 
return of 82% of patients to duty, resulting in very few 
evacuations (11 casualties per day) , thus reduceing the demand 
on combat transportation resources. An additional benefit is 
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that, by evacuating fewer casualties out of the combat zone, 
fewer replacements are required to be sent into the area. 

However, to operate under this policy, the medical 
battalion must be able to sustain an average number of 478 
occupied beds. This bed utilization could pose an immense 
burden on the supply system, besides making the medical 
battalion highly immobile. Moreover, this evacuation policy 
allows only a remote chance of accommodating a casualty spike 
of 120 patients. The consequences of this incapacity to 
handle incoming casualties are likely to outweigh the 
advantages of a high rate of return to duty. 

2. Assuring the Accommodation of Casualty Spike 

On the other hand, if the medical planner wants a 100% 
probability of accommodating the estimated casualty spike, he 
would choose an evacuation policy of 10 days. This evacuation 
policy results in an average of 335 occupied beds. This lower 
bed utilization results in a more mobile medical battalion 
requiring much less supply support to sustain its operation. 

Unfortunately, this evacuation policy diminishes the 
proportion of casualties returned to duty to 58%. Along with 
this reduction comes an increase in casualty evacuations to 26 
per day. This is double the number of daily evacuations under 
the 17-day evacuation policy. Additionally, twice the number 
of inexperienced reinforcements are called upon to enter the 
area. Each of these increased flows between the medical 
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battalion and higher echelons increases the demand on 
transportation resources. 

Clearly the medical planner has to decide which of 
these, and other, quantitative tradeoffs must be made for his 
specific scenario. The models developed in previous chapters, 
and illustrated by the example in this chapter, help the 
medical planner evaluate the evacuation policy tradeoffs. 



46 



VI. CONCLUSIONS 



This thesis has developed several analytical models by 
which the evacuation policy problem can be analyzed 
quantitatively. The study closes with some comments on the 
significance of the current model and suggestions for areas of 
further research. 

A. POTENTIAL MODEL USEFULNESS 

The analytical models developed in this thesis provide the 
medical planner with a basis for quantitatively evaluating the 
effects of various evacuation policies. These models can be 
used in either the planning or execution of a military 
operation. For planning purposes/ estimated parameters may be 
input into the model to evaluate medical planning strategies. 
Operationally/ the model can be updated dynamically with 
current estimates of the input parameters. This allows 
military commanders to evaluate the current policy and guides 
them in the selection of a new policy. 

The results for this study were calculated using both a 
hand-held calculator and the PC-based spreadsheet Excel. In 
military applications/ all calculations could be made on a 
programmable calculator/ programmed in any PC-based language, 
or implemented in most popular spreadsheets . 
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B. RECOMMENDATIONS FOR FURTHER STUDY 



The assumptions made in this study are reasonable for the 
Low-Intensity Conflict scenario. Complete evaluation of the 
significance of these assumptions is left to further research. 
There are also some specif ic ^ enhancements recommended to the 
models developed in this study. Model development assumed the 
same evacuation policy and evacuation schedule for all 
casualty types. The models could easily be augmented to allow 
different parameters for each casualty type. 

Casualty arrivals are modeled as Poisson for this study. 
The models could be changed to model group arrivals, thus 
incorporating extra variability in the Poisson arrival 
process . 

The model for the casualty treatment time distribution 
developed in this study is based on a single data set. A more 
robust model might be considered. A possibility is to develop 
an aggregate model based on a variety of specific disease or 
injury types. Additionally, the treatment distribution model 
could be generalized by studying various other casualty data. 

Another area for study concerns the modeling of the triage 
process. This study considered the decisions made by the 
triage officer to be perfect. Further research might apply 
decision theory to model imperfect triage, and thus study the 
uncertainty in the triage decision making. 

Finally, some enhancements could be made to account for 
inaccuracies in estimates of casualty rates, proportions and 
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treatment times . A general Bayesian model is suggested to 
allow planners to update these estimated parameters as actual 
data is obtained during the operation. 
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LOGNORMAL DENSITY FUNCTION, N=1827 LOGNORMAL PROBABILITY PLOT, N=1827 



APPENDIX A. AGSS DATA ANALYSIS PLOTS 
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Lognormal Fit for WIA Treatment Time 

LOGNORMAL DENSITY FUNCTION, N=358 LOGNORMAL PROBABILITY PLOT, N=358 
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APPENDIX B. PASCAL SOURCE CODE 



Program Evaluate_Erlang ; 

{ Author : Tracy L. Howard } 
{ Date : March 1993 } 
{ Use : Thesis Research } 
{ This program evaluates the Erlang-B distribution model } 
{ developed in this thesis for comparison with normal } 
{ approximations. } 



{ The mean treatment time data is calculated in Excel and} 
{ read from an input file. The output is analyzed in Excel.} 

{ The math unit is used for Factorial and Power functions.} 
Uses Math; 



Const 

MAX_BED_S I Z E = 1000; 
MAX_SCENARIO_DAYS = 90; 



Type 

DataArray = Array [ 0 . . 10*MAX_SCENARIO_DAYS] of Real; 



Scenario_Record = Record 
Evac 
NumBeds 
End; 

Cas_Record = Record 

Lambda , 

MST 

FreqAdmitted 

Data 

End; 



: Real; 
: Word; 



Real ; 

Word; 

DataArray; 



BedArray = Array [ 0 . .MAX_BED_SIZE] of Extended; 



Var 



Current 


: Scenario_Record 


CumDist , 
BedDistI 


: BedArray; 


WIA, 

Combined 


: Cas_Record; 


RTD, 




Mu 


: Real; 
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Counter, 

I 



: Word; 



DataOutFile, 
WIAInFile : Text; 



Procedure InitializeFiles (var WIA_Input, Output: Text); 

{ This procedure assigns the filenames for the input and} 
{ output files, then prepares the input file to be read } 
{ and output file to be written. } 

Var I : Word; 

Begin 

{ Assigns output file that will be written in comma } 

{ delimited format. } 

Assign (Output, 'd:\thesis\dat\output.csv'); 

Rewrite (Output) ; 

{ Assigns input file containing the means service times } 
Assign (WIA_Input, ' D : \THESIS\DAT\meansimm. txt ' ) ; 

Reset (WIAInFile) ; 

Readln (WIAInFile) ; 

End; {InitializeFiles} 



Procedure ReadData (var Data_File : Text; var Data : 

Data Array) ; 

{ This procedure reads the first 300 mean service times} 
{ from the input file. } 

Var Index : Word; 

Begin 

{ Initializes the data array to zero.} 

FillChar (Data, SizeOf (Data ) , 0); 

{ Reads the Data in the form of one item per line.} 

For Index := 1 to 300 do 

Readln (Data_File, Data [ Index] ) ; 

End; {ReadData} 



Procedure Erlang (NumberOf Beds : Word; 

Lambda, Mu : Real; var Prob : BedArray) ; 

{ This procedure calculates the probability distribution } 
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{ for the Erlang-B model. Note that this calculation } 

{ involves calculations of the form (X to the nth) divided } 

{ by n! . This algorithm evaluates these terms by brute } 

{ force (i.e., Very large numbers are needed). } 

{ Numerical methods could be used to develop a more } 

{ efficient algorithm. } 

Var I : Word; 

Rho, 

Numerator, 

Denominator : Extended; 

{ Extended type allows Real numbers with } 

{ Exponents in the range: -4932 to 4932. } 

Begin 

{ Initializes the probability array.} 

FillChar ( Prob, SizeOf (Prob) , 0 ) ; 

Rho := Lambda /Mu; 

Denominator : = 1.0; 

For I : = 1 to NumberOfBeds do Denominator := Denominator 

+ Power (Rho, I ) /Factorial ( I) ; 

For I : = 0 to NumberOfBeds do begin 

Numerator := Power (Rho , I ) /Factorial ( I ) ; 

Prob [I] := Numerator / Denominator; 

End; {For 1} 

End; {Erlang} 



Begin { Analytical_Model } 

{ Bed capacity of the facility} 

Current .NumBeds := 540; 

InitializeFiles (WIAInFile, DataOutFile) ; 

ReadData (WIAInFile, WIA.Data) ; 

{ Total casualty arrival rate} 

WI A. Lambda := 100 ; 

{ This loop evaluates the Erlang-B distribution for } 

{ evacuation policies from 0 to 30 days in 1/10 of a day} 
{ increments. } 
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For Counter ;= 1 to 300 do begin 
Current. Evac := Counter/ 10; 

Mu := 1/Wia . Data [Counter ] ; 

Erlang (Current . NumBeds , WIA . Lambda , Mu, BedDist I ) ; 
FillChar (CumDist, SizeOf (CumDist) , 0} ; 

CumDist[0] := BedDistI[0]; 

{ This loop writes the cumulative probability } 

{ distribution to the output file. } 

Write (DataOutFile, (WIA . Lambda/ (Mu*Current . NumBeds )) , ' 

CumDist [ 0] ) 

For I := 1 to Current .NumBeds do begin 

CumDist [I] := CumDist [1-1] + BedDistI[I]; 

Write (DataOutFile, ' , ' , CumDist [I] ) ; 

End; 

Writeln (DataOutFile) ; 

End; {For Current} 

Close (DataOutFile) ; 

Close (WIAInFile); 

End. {Evaluate_Erlang} 
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APPENDIX C. EXCEL WORKSHEETS 



The pages of this appendix contain printouts of 
Excel data worksheets used to calculate the results 
study . 



selected 
for this 
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Evac Casualty Arrival Lognormal Parameters 

P(Admission) Sched Rate X ay Mean Std Dev 

DNBI 0.94974093 4 3^5 52l 0.737 2A 3.74935753 9.100604442 

WIA 10.59777228 4 | i L5 22.5 0.751 2.25\ 4.0806243351 10,92187251 

Total 75 
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E[S] Total E[N] Normal Appr 0.95 P RTD 

3.371507 252.863 260.9810887 0.002012 

3.330691 249.8018 264.2010754 0.026187 

3.283696 246.2772 267.9097289 0.081858 

3.283696 246.2772 267.9097289 0.157569 
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